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Recent work has indicated that ghost imaging may have applications in standoff sensing. However, 
most theoretical work has addressed transmission-based ghost imaging. To be a viable remote- 
sensing system, the ghost imager needs to image rough-surfaced targets in reflection through long, 
turbulent optical paths. We develop, within a Gaussian-state framework, expressions for the spatial 
resolution, image contrast, and signal-to-noise ratio of such a system. We consider rough-surfaced 
targets that create fully developed speckle in their returns, and Kolmogorov-spectrum turbulence 
that is uniformly distributed along all propagation paths. We address both classical and nonclassical 
optical sources, as well as a computational ghost imager. 



I. INTRODUCTION 

Ghost imaging is a transverse active-imaging technique 
that exploits the correlation between two light beams to 
image a target without spatially-resolving measurements 
of the light beam that has undergone target interaction. 
The two beams are the signal, which interacts with the 
target and then is measured by a single-pixel bucket de- 
tector, and the reference, which is directly measured by 
a high spatial-resolution detector. An intensity cross- 
correlation between the signal beam that encounters the 
target and the reference beam that impinges on the high 
spatial-resolution detector imparts target information — 
the ghost image — to the cross correlation between the 
photocurrents obtained from the two detectors, informa- 
tion that is unavailable from either photocurrent alone. 

Most ghost-imaging experiments and related theory 
deal with the transmissive case, in which the bucket de- 
tector is placed behind the target, and we image the 
target's intensity-transmission profile. In the first ghost 
imaging experiment, the signal and reference beams 
were the quantum-mechanically entangled outputs from 
a spontaneous parametric downconverter (SPDC) pQ, 
which afforded them a phase-sensitive cross correlation 
stronger than permitted by classical physics [2]. The 
system was run at low flux, in which the SPDC's post- 
selected output state, within the photodetectors' re- 
sponse time, could be taken to be a biphoton; the ghost 
image was then formed by counting coincidences between 
the bucket and high spatial-resolution detections instead 
of a photocurrent cross correlation. Subsequently, ghost 
imaging was done with a high- flux pseudothermal source, 
i.e., a classical-state source realized by sending a laser 
beam through a rotating ground-glass diffuser and a SO- 
SO beam splitter to create identical, spatially-incoherent 
signal and reference beams with a phase-insensitive cross 
correlation. The ghost image was then obtained by 
continuous-time cross correlation of the bucket and high 
spatial-resolution photocurrents [31 |3]. Later, the the- 
ories of SPDC and pseudothermal ghost imaging were 
united in a Gaussian-state treatment that showed, in 
both cases, that ghost-image formation by photocurrent 
correlation arises from classical coherence propagation. 



but the nonclassical SPDC source offers much higher im- 
age contrast and a modest spatial-resolution advantage 
in near- field operation [21 [Sj [6] . 

Recently, a computational ghost-imaging method was 
introduced, in which the need for a reference beam is re- 
moved [7J IB] • Conventional ghost imaging relies on the 
intensity cross-correlation between the signal beam on 
the target and the reference beam on the high spatial- 
resolution detector. In SPDC and pseudothermal ghost 
imaging, these beam patterns are random. If, how- 
ever, a spatial light modulator is used to impart a se- 
quence of known spatial patterns to a laser beam prior 
to propagation to the target, then diffraction theory 
can be used to calculate noiseless versions of the asso- 
ciated reference-beam intensity patterns. These com- 
puted reference-beam results can be employed, for cross- 
correlation ghost-image formation, as if they came from 
photodetection of a physical reference beam. More im- 
portantly, computational ghost imaging opens the door 
for more advanced computational image reconstruction 
techniques, notably compressive sensing [Sj. 

All of the preceding discussion has been framed for 
ghost imaging of transmissive targets. Recent experi- 
ments have shown the feasibility of an alternate ghost- 
imaging configuration, in which the bucket detector views 
the target in reflection, rather than in transmission |10) . 
In this reflective setup, the source and detector can, in 
principle, be moved arbitrarily far away from the target, 
suggesting that ghost imaging could be a viable remote- 
sensing technique. Atmospheric turbulence will surely 
be a relevant concern in any such standoff-sensing appli- 
cation, just as it is for astronomical imaging and laser 
radar. To date, there has been theoretical study of the 
impact of turbulence on ghost imaging in transmission 
[11) . but there has only been an initial turbulence-free 
theory development for ghost imaging in reflection [12| . 
and that work was limited to pseudothermal sources. 

In this paper we will extend the analysis from [12j 
to include SPDC ghost imaging in reflection and, for 
both SPDC and pseudothermal operation, the presence 
of atmospheric turbulence in the propagation paths. The 
setup we shall consider is the lensless ghost-imaging con- 
flguration shown schematically in Fig. [l] [13]. Pseu- 
dothermal ghost imaging can be described quantitatively 
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using semiclassical photodetection theory — in which the 
fields are described classically and photodetection incurs 
the shot noise associated with the discreteness of the elec- 
tron charge — whereas quantum photodetection theory is 
required for the SPDC case, owing to the nonclassical 
nature of its output state [2]. Nevertheless, to provide a 
unified development, we will employ quantum photode- 
tection to characterize both pseudothermal and SPDC 
ghost imaging. We begin, in Sec. [11} with the theoretical 
framework for our analysis. Here we will describe our 
models for the pseudothermal and SPDC sources, reflec- 
tion from the target, atmospheric propagation, quantum 
photodetection, and image formation. In Sec. Ill we will 
derive expressions for the spatial resolution and image 
contrast for pseudothermal and SPDC ghost images, and 
in Sec. |IV| we will derive their respective signal-to-noise 
ratios (SNRs). Then, in Sec.|v) we shall consider a com- 
putational ghost imager, in which the reference beam is 
calculated instead of measured [7]. Because the refer- 
ence beam is not an optical field, we find it convenient in 
Sec. |V] to employ the semiclassical theory. This section 
will also include a quantitative performance comparison 
between pseudothermal, SDPC, and computational ghost 
images formed in reflection. We conclude, in Sec. VI 
with a discussion of our findings, and an assessment of the 
situations in which reflective ghost-imaging might have 
practical value. 
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FIG. 1: (Color online) Light from a spatially-incoherent op- 
tical source — either a type-II phase matched SPDC or a laser 
rendered spatially incoherent by passage through a rotating 
ground-glass — is separated into signal and reference beams 
by a beam splitter. For the pseudothermal (laser) case, this 
is a 50-50 beam splitter. For the SPDC, this is a polarizing 
beam splitter. The reference beam travels L m to the CCD 
camera, while the signal beam travels Lm to the target, and 
the reflected light travels Lm to the bucket detector. The 
image is then formed by cross correlation between the pho- 
tocurrents obtained from the bucket detector and the high 
spatial-resolution (CCD camera) detector. 



II. SETUP 

A. Optical Sources 

We are interested in ghost imaging with classical 
pseudothermal light, and the nonclassical output of 
an SPDC. In both cases our signal and reference 
beams will be represented quantum-mcchanically as 
y/ photons / m^s-units , -j-z-going, positive-frequency, field 
operators — Es{p,t)e~^'^°* and £^h(p, t)e~*"°*, as func- 
tions of transverse-coordinate vector p and time — at the 
output of the Fig. [T] beam splitter. Their baseband field 
operators have the canonical commutation relations 

[Ejip,,h),EK{P2,t2)]=0, (1) 

[Ejip,,h),El,{p2,h)]^ 5jkS{p,- p2)S{t, -h), (2) 

for J = S,R and K = S,R, where Sjk is the Kro- 
necker delta, and S{-) is the unit impulse. Following [5], 
we will take both the pseudothermal and SPDC sources 
to produce signal and reference beams that are in zero- 
mean, jointly Gaussian states, which are then completely 
characterized by their nonzero autocorrelation and cross- 
correlation functions. All the autocorrelations will be as- 



sumed to have the same Gaussian-Schell model form, 

{El{puti)EK{p2,t2)) 

2P |Pll^ + |P2l^ PI-P2I'' (tl-*2)^ 

= — s-e "5 e ^''0 e ^^0 , (3) 

for K = S, R, where P is the photon flux, ag is the 
source's intensity radius, po is its coherence length, 
and Tq is its coherence time. The distinction between the 
Gaussian-state models for the pseudothermal and SPDC 
sources lies in their signal-reference cross correlations, as 
we will now explain. 

The pseudothermal signal and reference beams start as 
continuous- wave (cw) laser light, that is passed through 
a rotating ground-glass diffuser to make it spatially in- 
coherent. Signal and reference beams are then obtained 
by 50-50 beam splitting. The signal and reference thus 
have no phase-sensitive cross correlation, and the maxi- 
mum phase-insensitive cross correlation allowed by clas- 
sical physics given their autocorrelation functions, i.e., 

{Es{pi,ti)ERip2,t2))c ^0. (4) 

{El{puti)ER{p2,t2))c 

op IPll^ + IP2l^ Pl-P2l^ (tl-*2)^ (5) 

= — =-e e ■''•0 e , 
Tra^ 
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where the subscript C indicates that these are classical- 
state cross correlations. The coherence time is directly 
related to how fast the ground-glass rotates, the coher- 
ence length is typically on the order of a few wavelengths, 
and the intensity radius is set by the radius of the initial 
laser beam. 

A type-II phase-matched, collinear, cw SPDC source 
emits orthogonally-polarized, co-propagating signal and 



idler beams, which are separated by a polarizing beam 
splitter to become the signal and reference beams for 
ghost imaging. They are in a maximally-entangled, 
jointly-Gaussian state with no phase-insensitive cross 
correlation, and the maximum phase-sensitive cross cor- 
relation permitted by quantum theory given their auto- 
correlation functions, namely [5] 



{Es{piM)Er{p2M))c 



2P IPlI +|P2 

— s-e "0 



PI -P2 I 

2^? 



_ (tl-t2)" 




(tl-*2)" 



(6) 
(7) 



where the subscript Q indicates that these are quantum 
(nonclassical) cross correlations. 



B. Propagation through Turbulence 

The signal and reference beams in Fig. [l] first propa- 
gate L m over spatially-separated paths from their shared 
source location [13^. The reference beam is then mea- 
sured by a high-spatial resolution CCD camera, while 
the signal field reflects off the target and travels Lm 
to the bucket detector. We will assume statistically- 
independent turbulence on all three paths: the reference 
path (R) from the source to the CCD; the signal path 
{S) from the source to the target; and the target-return 
path (T) from the target to the bucket detector. 

To account for the Kolmogorov-spectrum turbulence 
that is uniformly distributed along the three propaga- 
tion paths, we make use of the quantum version of the 
extended Huygens-Fresnel principle [HI [T3] 

E'^{p\t) = jdpEUp^tf-^^^—^ e^'"('''''), 

(8) 

where /cq — ^^o/c ~ 27r/Ao, and tp^ip' i p) is a complex- 
valued random process that encapsulates the effects of 
turbulence on path m from p to p' . The real and imag- 
inary parts of tpmip', p), which we will denote Xm{p' , p) 
and (f),n{p',p), are the logamplitude and phase fluctua- 
tions imposed on the field received at transverse coordi- 
nate p' over path m from a point source at transverse 
coordinate p. Because all the path lengths in Fig. [l] are 
equal, we have suppressed the L/c line-of-sight propaga- 
tion delay in Eq. We have also suppressed any time- 
dependence of the turbulence. Physically, this amounts 
to assuming that the integration time, T/, employed in 
the Fig. [T] correlator is shorter than the ~ms coherence 
time of the turbulence. However, that assumption is not 



necessary for deriving the spatial resolution and image 
contrast achieved in the presence of turbulence, although 
it will be relevant to our signal-to-noise ratio analysis. 
The mutual coherence function of e'^™''' -^^ is taken to 

be [niin] 

^ g-(|p'l-P2l^ + (pi-P2)-(Pi-P2) + |pi-p2|^)/2pJ„^ ^g^ 

which is the square-law approximation to the rigorous 
5/3-law behavior for uniformly-distributed turbulence 
[18] . In Eq. is the turbulence coherence length 

on path m, given by 



(10) 



in terms of the constant turbulence strength-parameter, 
,„, along path m [16) . 
We will assume that the path length L satisfies the far- 
field coherence-propagation conditions for pseudothermal 
and SPDC light discussed in ^,6j, viz., fcoaoPo/2i ^ 1, 
and koal/2L <^ 1, respectively. For convenience, we also 
define pL — iL/k^aQ and = 2L/kQpo. These are the 
far-field, turbulence-free, on-target coherence length and 
intensity radius produced by the Gaussian-Schell auto- 
correlation function from Eq. ([s]). 



C. Target Reflection 

Most real-world targets have surfaces that are suffi- 
ciently rough — on the scale of an optical wavelength — 
that reflections from them present fully-developed 
speckle patterns at appreciable distances from the sur- 
face [19]. We will model such a target as a planar 2- 
D object with random, microscopic height variations — 
having depths on the order of several wavelengths, and 
transverse variations on the order a wavelength — and an 
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average intensity-reflection coefScient T{p). Thus, fol- 
lowing laser-radar theory |20J, we shall assume that this 
target has a random field-reflection coefScient, T{p), that 
is a zero-mean, complex-valued Gaussian random process 
characterized by the autocorrelation function 



{T*{p,)T{p2)) = \lT{pi)d{pi ~ p2), 



(11) 



where the deterministic pattern T{p) is what we are tying 
to image. 

Target reflection is a passive process, so we require 
|2^(p)| < 1, which is in conflict with the Gaussian statis- 
tics and the delta- function term in Eq. (111. That delta 



function, however, leads to quasi-Lambertian reflection, 
implying that at standoff distances a realistic bucket de- 
tector will only capture a very small fraction of the re- 
flected light. Hence using the reflection statistics from 
the preceding paragraph is a reasonable approximation 
for ghost imaging of rough-surfaced targets. However, 
because some light may be absorbed by or transmitted 
through the target, the quantum model for target reflec- 
tion needs to be 



ET{p,t)=T{p)&sip,t) 



^l-\T{p)\^E,ac.s{p,t) 

(12) 



Here E^ac,s{Pi t) is a vacuum-state field operator, needed 
to preserve the free-field commutator relations given, for 
the source fields, in Eqs. ([T]) and ([2|. 

In what follows we will assume that the target's 
intensity-reflection coefficient, T(p), is space-limited 
to the region illuminated by the signal beam, i.e., 
|T(p)| = for IpI > a^. This condition will give the 
ensemble-average pseudothermal and SDPC ghost im- 
ages convolution-integral relationships with T(p) from 
which to determine their spatial resolutions. See [H |5] 
for information about ghost imaging field-of-view when 
the preceding restriction on T(p) is not satisfied. 



CCD and the bucket detector, formalized as 

1 r'^i/'^ 



G{Pp) 



Tj 



Ti/2 



dt ip{t)ib{t), 



(14) 



where Pp is the center coordinate of the CCD's pth pixel, 
and Ti is the correlation integration time. From quantum 
photodetection theory [22], we have that 



i„i{t)=q dTPmiT)h{t-T), m=p,b 



where 



Pm{t) = / dpEUp,t)Era{p,t) 



(15) 



(16) 



is the effective (sub-unity quantum efficiency) photon flux 
on the active region Am for detector m, and h(t) is the 
impulse response of the detectors' output circuit, which 
we have taken to be the same for both the bucket detector 
and the CCD camera. 

To enable the performance evaluations in the sections 
to follow, we make the following assumptions. First, each 
pixel is small enough that the the average photon flux 
impinging on it obeys {Pp{t)) = Ap{El{pp,t)Ep{pp,t)) , 
where Ap is the area of the Ap. Second, we assume that 
the detector has a Gaussian impulse response. 



hit) 



n 



Stt 



(17) 



whose bandwidth, Qbi greatly exceeds that of the in- 
cident hght for the pseudothermal source {^bTq » 1), 
and is much smaller than that of the incident light for the 
SPDC source {ilsTo <C 1). Finally, we assume that the 
correlation integration time is sufficiently long to capture 
a ghost image, so Tj ^ Tq for the pseudothermal source, 
and Ti > rj^i for the SPDC source. 



D. Photodetection and Image Formation 



III. RESOLUTION AND CONTRAST 



Loss also occurs in the photodetection process, be- 
cause the detectors in Fig. [T] will, in general, have sub- 
unity quantum efficiencies. Therefore, the field operator 
Ep{p,t) driving the photocurrent, ip(t), from CCD pixel 
p, and the field operator Ef){p,t) driving the photocur- 
rent, ib{t), from the bucket detector satisfy pH [22] 



Em{p,t) ^ ^E[{p,t) + ^l^Kac/{p,t), (13) 

for {m,(.) = (p, i?) or (to,£) = {b,T). Here, r] is the 
quantum efficiency, which is taken to be the same for both 
the bucket detector and the CCD camera, and Eyac,e{p, t) 
is a vacuum-state field operator whose inclusion models 
the loss incurred by having 77 < 1, cf. Eq. (12). 



We will assume ghost-image formation via continuous- 
time, pixel- wise correlation of the photocurrents from the 



The spatial resolution and contrast of the ghost image 
are properties of the average cross correlation, {C{pp)). 
To determine this ensemble-average cross correlation we: 
(1) write ib{t) and ip{t) in terms of the field operators 
that drive them, Eh{p,t) and Ep{p,t); (2) apply the 
canonical commutator relations to normally-order the 
resulting photocurrent-product expressions; and (3) use 



Eq. (13) to write the detected fields in terms of E'j^{p,t) 
and E'rp{p,t). From the linearity of expectation, we can 
bring the ensemble-average operation inside the time- 
domain integration to get 



{C{pp)) = q^r/Mp j dn J dT2 dph{t - Ti)h{t - T2) 

X {&l{pp,n)&^{p,T2)E'R{Pp,Ti)E'T{p,T2)). (18) 
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We then use Eqs. ([8]) and ( 12 ) to back-propagate the are 
target-return field operator to the target, and find 

{E'lip^, n)E'^ip, T2)E'j,{pp, n)E'Ap, T2)) 

- Jdp2j dp3 {T*{p2)T{p^)){e^T{P,P2)^^Ti.P.P.)^ 

X (pp , n ) i^s^ (p2 , T2 ) i;^ (Pp , ri ) i;^ (P3 , T2 )) 

p-iko(\p-2f ~\p3f) /2L+ikap-(p2- Ps) / L 



(19) 



where we have used the statistical independence of the 
target's reflection coefficient, the signal and reference 
fields at the source, and the turbulence on all three prop- 
agation paths. 



Substituting for (T* (^2)^(^3)) from Eq. ( 11 ) and per- 



forming the P2 integral eliminates the turbulence term 
from the preceding (C(pp)) expression, which shows that 
turbulence on the the target-to-bucket path has no im- 
pact on the ghost image's spatial resolution or image con- 
trast. Now, by back-propagating the field operators to 
the source — by means of Eq. (|8| with koaQ/2L <^ 1 — we 
have that 

{E'^ (Pp , n )E'^ (p, T2)E'j, [pp ,n)E'j.{p,T2)) 

= X^ldp2 T{P2) I dp[ I dp'2 I dp'i I dp'^ 

X {El{p[,T,)El{p'2,T2)ER{pln)Es{P2,r2)) 

V /pi'RiPp^P'l) p4'RiPp,Pl)\ /p4's(P2^P'2) pi's(P2:P'2)\ 



X e 



iko[pp-(p[-Pi)+P2-{p'2-p'2)]/L 



(20) 



where we have again invoked the statistical independence 
of the signal and reference fields at the source and the 
turbulence on the source-to-CCD, and source-to-target 
paths. 

We can now directly evaluate all remaining moments. 
The turbulence moments are available from Eq. (|9]), 
and the fourth-order field moment can be expressed in 
terms of second-order moments with the Gaussian-state 
moment-factoring theorem |23j : 

{El{p[,n)El{p'2,T2)ER{p'l, Tl)Es{P2, r2)) = 

(4(pi, ri)4(p^, T2)) {EM, T,)EsiP2, r2)) + 
(4(p'i, ri)^«(p'/, ri))(4(p^, r2)^s(p'2', ^2)) + 
(4(p'i, ri)^s(p^', T2)) {eI{p'2,T2)Er{pI n)). (21) 

Until this point, our analysis applies equally well to the 
pseudothermal and SPDC ghost imagers. That equiva- 
lence disappears, however, when we evaluate the second- 
order moments in Eq. (21 1 using Eqs. ([3|-([5|) for the 
pseudothermal ghost imager, and using Eqs. ([sj^ ([6]), and 
^ for the SPDC ghost imager. The results we obtain 



{C{Pp))c = 



L2 



dpr(p) 



^-\p-Pp\ /aPL 



(22) 



for the pseudothermal (classical) imager, and 



(C(Pp))q = 



L2 



dpTi-p) 



Q,bTq e 



-Ip-Ppl'/aPl 



1 



(23) 



for the SPDC (quantum) ghost imager, where A^, is the 
area oi Ah, 1 = PToPq/aq — PTop\/a\ is the source 
brightness in photons per spatiotemporal mode, and 



(24) 



is a turbulence-induced performance-degradation factor. 
We see that both average correlation functions are the 
sum of the same featureless background term 



q^rf-ApAy, 
L2 



2P 



dpT{p), 



(25) 



and different image-bearing terms 
aL2 i 



Ci,c{Pp) 



dpT{p)e 



\P-Pp\ 



(26) 



for the classical-state (pseudothermal) source, and 



Ci,q(Pp) 



q'r^^ApAb 



ai2 



2P 

Trn^ 



dpT{-p)e 



\p-pp\' 



(27) 



for the nonclassical-state (quantum) source. 

The average pseudothermal ghost image is erect, and 
consists of a scaled version of the target's intensity- 
refiection coefficient T(p) convolved with a Gaussian 
point-spread function whose e^^ radius is pL\/ot. The 
average SPDC ghost image is inverted, with a differ- 
ent scaling but the same point-spread function. Thus, 
both the pseudothermal and SPDC ghost images have 
the same spatial resolution, 



XoL 



XqL /. a? 



H-^(1.09ft2L)f (C2_^)f +(C2^^)^ 



(28) 
.(29) 
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where we have utihzed Eq. (10 1 to write the turbu- 



lence coherence lengths on each path in terms of their 
respective strength parameters. Equation (28) shows 



that turbulence-limited resolution prevails for the far- 
field lensless ghost image when either the source-to-target 
coherence length (ps) or the source-to-reference coher- 
ence length {pr) becomes smaller than the source size 
(ao). When ps ^ uq and p/j ^ ag, far-field lensless 
ghost imaging in reflection yields the same XoL/irao spa- 
tial resolution previously reported for far-field lensless 
ghost imaging in transmission [51 [T^]. However, even 
when turbulence does not degrade ghost-image spatial 
resolution, there is still a significant difference between 
the transmissive and reflective cases, namely the Ab/L'^ 
factor that is present in our results but is absent from 
those in [5]. This factor, which will obey A^/L'^ <C 1 in 
a standoff-sensing scenario, is the bucket detector's an- 
gular subtense at the target, and represents the fraction 
of the quasi-Lambertian reflected light that this detector 
collects. No such factor appeared in [5], because that 
work assumed the bucket detector collected all the light 
that was transmitted through the target. 

The featureless background terms in our pseudother- 
mal and SPDC ghost images affect the images' contrast. 
To assess these contrasts, we will adopt the definition and 
approach presented in In particular, the image con- 
trast C is the ratio of the difference between the brightest 
and darkest pixels in the image to the featureless back- 
ground: 



C 



maxT^ [C'(pp)] - miuTzjC (Pp 
Co 



(30) 



where TZ is the region within which T(p) is nonzero. We 
will assume that the ghost imager completely resolves the 
target, so that we can use the approximation 



dpT{±p)e 



\p-pp\' 



■nap\T{±p). 



(31) 



Defining At — JdpT{p) as the target's reflective cross- 
section, and taking max-jj [T(p)] = 1, we get the following 
image contrast expressions for our two ghost imagers, 



At ' 



and 



Co = 



At 4 



1 



4^/^2: 



(32) 



(33) 



If |T(p)| ~ 1 where it is nonzero, then At is the area 
of the target. Moreover, p\ is the spatial-resolution area 
in the absence of turbulence. It follows that the pseu- 
dothermal ghost-image contrast is approximately l/iV, 
where N is number of no-turbulence resolution cells in 
the image. Note that our analysis to this point has as- 
sumed a DC-coupled correlator in Fig. [l] The contrast of 



a pseudothermal ghost image can be improved by sup- 
pressing the background term through AC-coupling, i.e., 
measuring the cross covariance, rather than the cross cor- 
relation, between the bucket and high spatial-resolution 
photodetector's outputs, as was done in the pseudother- 
mal experiments of Scarcelli et al. [24] . 

The pseudothermal imager's contrast does not depend 
on its source brightness, I, but that for the SPDC does. 
At high source brightness, 1^1, the SPDC image con- 
trast reduces to the result we found for the pseudother- 
mal case. The normal operating regime for a cw SPDC, 
however, is one of low brightness, 1^1, for which 
Eq. ( 33 ) reduces to 



Cc 



(34) 



The number of photon pairs emitted by an SPDC in one 
detector time-constant is low, so fts/P ^ 1, and our 
assumption that the target is contained within the region 
illuminated by the signal beam implies that aj^/AT > 1 
if |T(p)| ~ 1 where it is nonzero. Consequently, unlike 
the pseudothermal case, Cq ^ 1 is the norm for SDPC 
ghost imaging in DC-coupled operation. 



IV. SIGNAL-TO-NOISE RATIO 

In keeping with the analysis done in [5] , we will be eval- 
uating the signal-to-noise ratio (SNR) pixel-wise as the 
ratio of the square of the expected value to the variance 
for the photocurrent cross-correlation function. 



SNR: 



(CiPp))' 



{C{Pp)? 

Var[C(p,)] {C^pp)) - {C{pp)Y 



(35) 



The main complication in evaluating Eq. ( 35 ) is the sec- 
ond moment. 



Ti/2 Ti/2 



(C'iPp)) 



dt 



du {ip{t)ib{t)ip{u)ii,{u)) . (36) 



-Tj/2 -Tj/2 



After using Eq. ([l5| to write the currents in terms of 
the fields, Eqs. (|l| and ^ to normally order the fields, 
and Eqs. ([s]) and ( [l2| ) to propagate the fields back to 
the source, we are left needing to evaluate an eighth- 
order moment of the fields, two sixth-order moments of 
the fields, a fourth-order moment of the fields, a fourth- 
order moment of the target, three fourth-order moments 
of the turbulence, and propagation of these terms over 12 
unique paths. To make the analysis more tractable, we 
use several simplifying assumptions. First, we AC-couple 
both of the detector outputs by inserting a narrowband 
notch centered at the origin in frequency space, resulting 
in the frequency response 



(37) 
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where i^N ^ ^b- We also assume that the ghost imager 
can resolve all structure on the target, allowing us to use 
the approximation in Eq. (31). 

To evaluate the fourth-order turbulence moments, we 
can no longer use the square-law approximation in Eq. ([9| 
that models the turbulence's mutual-coherence function 
as Gaussian, because this approximation is not valid for 
fourth-order moments j25j . To evaluate these terms we 
assume that the turbulence coherence length on each 
path is much larger than the on-target coherence length 
PL- We also assume the turbulence coherence lengths 
are all much larger than oq, implying that the turbu- 
lence functions are independent of coordinates at the 
source and detection planes. Although these assump- 
tions are somewhat limiting, they still allow weak-to- 
medium strength turbulence for a typical situation. For 
instance, if the target is standing off 1 km, the source 
aperture has radius oq = 3 cm, and we are at a wave- 
length Xq — 1.5 /im, this assumption is satisfied for 



< 10- 



,-2/3 



Moreover, this situation has no 



spatial-averaging of the turbulence at the target, and 
thus represents a worst-case scenario for the SNR. 

For p'l and being coordinates at the CCD or bucket 
detectors, and pi being coordinates at the target, the 
preceding coherence-length assumptions imply that 



(^g'^'n (p'l .Pi ) g'Am (p'l' .Pi ) gV-i^i (P2 .P2 ) gV"™ (Pa .P2 ) ^ 

~ (g'/',*^ (o.Pi ) g'A™ (o.Pi ) g'A™ (o,p2) gi>™ (o,p2 ) \ 



„2Rc[V'™(0,pi)]+2Rc[iA„(0,p2)]\ 



(38) 



Decomposing the turbulence function into its log- 
amplitude and phase components as ipmipi, P2) ~ 
Xm{Pi,P2) + i(t>m{Pi, P2), and taking the logamplitude 
to be Gaussian distributed with mean equal to minus its 
variance [U], then gives 



^gV'm (p'l .Pl ) gV-m (p'l .Pi ) ^tp'm (p'i'P-i) g-iAm (P2' .P2 ) ^ 

~ p4i<:x™(pi-P2) 



(39) 



where 



^x™(Pi -P2) = (Ax™(0,pi)Ax™(0,p2)), (40) 

for = Xm ~ {Xm)i IS the log-amplitude covariance 

function. The logamplitude variance will be taken to be 
[TTj 

= K^^O) = 0.124 Cl^X^'L^'/^ (41) 



where we have used the Rytov-approximation expression, 
whose validity is ensured by our assumption of weak-to- 
medium strength turbulence. 

Even with these simplifications, the SNR evaluation is 
quite tedious, if somewhat straightforward. The higher- 
order field moments and the fourth-order target moment 
can be reduced to a sum of the products of second- 
order moments by application of the Gaussian moment- 
factoring theorem. The AC-coupling removes several of 
these terms, and we are left with integrals and Fourier 
transforms of the product of Gaussian terms, the tur- 
bulence correlation functions, and the target's intensity- 
reflection coefhcient, T{p). To simplify our final SNR 
expressions, we define two new terms, 

A't = j dpT^ip), (42) 

^=(4^/^^" e-HV20(«.,4y^), (43) 

where 0{C,D) is the two-circle overlap function for cir- 
cles of diameter D, 



0{C,D) = 



2 
0, 



COS 



-'(») 



D 



1 - 



D2 



for Id < D 
elsewhere, 



(44) 



/3 = Ab/naQ, and we have utilized the dimensionless dif- 
ference coordinates v = pLko{p' — p")/L, where p' and 
p" are coordinates at the bucket detector. 

Both A'rp and F have significant physical interpreta- 
tions. Similar to At, A'j, is another measure of the re- 
flective area of the target, and thus directly related to the 
number of on-target resolution cells, and the subsequent 
amount of time it takes to form an image. F is a met- 
ric for the spatial averaging of the target-induced speckle 
over the bucket detector. For very small detectors F w 1, 
while for large detectors it is inversely proportional to 
the area of the bucket detector, viz. F ^ 1/2/3. To fur- 
ther simplify our SNR expressions, we make two final 
assumptions: first that ^ 30, which is equivalent 

to saying that the ghost image consists of at least 10 x 10 
resolution cells; and second, that the bucket detector is 
at least as large as the area of the source beam, or /3 > 1. 
Applying these conditions we arrive at the following SNR 
expressions for pseudothermal and SPDC reflective ghost 
imaging through atmospheric turbulence. 
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SNRf 



To 
Ti 



A" Source + A^Path 



^bTq 2t 



Tf 



(45) 



-A^Detect + — A^Mix 
Ti Ti 



SNRq = 



1 



^bTi 



A^ Source + A^Path 



1 



2 W2^2 
Ti 



4.V^T 
A^Detect 



1 + 



To A 2 



Ti 



A^Mix 



2 1 



(46) 



where, for ease of comparison with the pseudothermal 
imager, we have inverted the image coordinates in the 
SPDC case. 

The terms that appear in the noise denominators of 
Eqs. (45) and (46) are as foUows: 



A Source = 



yl^(l + /3-i)e4('^H+'^l+4) 



A^Path = r^(pp) Te'^f'^H+'^l+'^T) (r + 1) - 1 

A^Detect = V ir, 
16^/2Ap772I2 Ab 

A^Mix^^^e^^?^ 



+ 



VlT'(Pp) 

Ap-ql 



„4(a|+<T|,) 



(47) 
(48) 
(49) 

(50) 



where for £ — R,S,T denotes the turbulence-induced 
logamphtude variance on the reference (R), signal (S*), 
and target (T) paths. Each of these noise terms has an 
important physical interpretation, as we explain below. 

The first noise term, A^Source, captures the fluctu- 
ations associated with decorrelation of transverse loca- 
tions in the image arising from the source's time evolu- 
tion. Thus, it is inversely proportional to the number of 
on-target resolution cells, p\/A'j,, and the number of dif- 
ferent source patterns measured. For the pseudothermal 
imager, for which ^bTq 3> 1 holds, this is the number of 
field-coherence times in the integration time, so we find 
A^Source scaled by Tq/T/. For the SPDC imager we 
have risTo ^ 1, so it is the detectors' response time that 
limits the measurements' correlation behavior, hence we 
find A^Source scaled by A/Q,bTi in that imager's SNR 
denominator. 

The A^Path term is due to the optical-path interac- 
tions that degrade the measurement, i.e., the turbulent 
atmosphere, and the scattering off the rough target. The 
noise contribution from target scattering can be miti- 
gated by aperture averaging, which, as noted above, is 
quantified by P. The effect of the turbulence is some- 
what more complicated, as it contributes a scaling factor 
e^'^i to all field-variation terms on path I. In addition, 
the logamphtude fluctuations also cause a deviation that 



results in the squared-mean term in the denominator of 
Eq. ( 35 ) no longer being canceled. These effects are par- 



ticularly devastating because A^Path is the only noise 
term whose strength is independent of the detector in- 
tegration time. This is because we have assumed time- 
independent target behavior and T/ less than the ^ms 
turbulence coherence time. Consequently, regardless of 
the source brightness and integration time, A^Path sets 
an ultimate upper limit on the SNR. 

The noise term A^Detect results from the statistics of 
photodetection at both the reference and bucket detec- 
tors, and as such is inversely proportional to the photon 
fluxes impinging on their active regions. Thus, A^Detect 
is inversely proportional to the size of each pixel, Ap, and 
the solid-angle subtense factor Ab/L'^. Detection noise 
is mitigated by temporal averaging, but exacerbated by 
having a fast detector, so the pseudothermal A^Detect 
is scaled by ^IbTq/Ti, while for the SPDC imager it is 
scaled only by Tq/Ti. The SPDC case is also scaled by 
by 1 -f l/4-y/7rI, a result of that source's stronger-than 
classical correlation. 

The A^Mix term contains noise contributed by mixing 
between source plus path noise on the reference arm with 
detection noise on the signal arm, and vice versa. For 
both imagers this term has a temporal scaling of Tq/Ti, 
and for the SPDC imager the quantum correlation im- 
parts an additional factor of 2/^/3 + l/Ay/nl. 

To get a better understanding of the SNR behavior for 
pseudothermal and SPDC ghost imagers, we will consider 
three limiting cases. The first is the long integration-time 
limit in which A^Path dominates the noise denominators 
for both the pseudothermal and SPDC ghost imagers, 
leading to the identical maximum (saturation) SNR value 



SNR, 



-4(<T^+<7|+<T^) 



(51) 



regardless of the strength of the illumination. 

For our next two cases, we will assume the integration 
time is sufficiently short that SNR saturation is not ap- 
proached. Then, in the high-brightness (I ^ 1) case, for 
which detection noise can be neglected, performance is 
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dominated by source fluctuations, A^Source, yielding 



Tj y2^p2 g-4K+<T|+^|) 



8 + 



T'iPp), (52) 



T'iPp), (53) 



for pseudothermal (C) and SPDC (Q) imagers, respec- 
tively. With all parameters other than source coherence- 
time being equal, the SPDC high-brightness SNR exceeds 
that of the pseudothermal source by a factor ^bTq » 
1, where Tq is the pseudothermal source's coherence 
time. This is because the much faster fluctuations in 
the SPDC output fields lead to faster decorrelation, and 
hence higher SNR. That said, however, the reader is re- 
minded that cw SPDC operation is ordinarily in the low- 
brightness regime, hence we turn our attention now to 
that special case. 

In the low-brightness, X ^ 1, regime, SNR behavior 
is dominated by the beat between the detection noises, 
A^Detect. Here we find that 



SNR.,Q = ^^r(Pp)^, 



(54) 
(55) 



for the pseudothermal and SPDC cases. The SPDC SNR 
is linear in source brightness, whereas the pseudothermal 
SNR is quadratic in that parameter. This behavior is 
to be expected, because the SPDC low-brightness out- 
put approaches the biphoton state, which is known to 
have lower detection noise than classical-state light. Fur- 
thermore, Tq for the pseudothermal source will greatly 
exceed Tq for the SPDC, and, as noted in the preced- 
ing paragraph, VLbTq 3> 1 holds for the pseudothermal 
source. It follows that the low-brightness SPDC SNR 
is much higher than that of the pseudothermal imager. 
The reader is cautioned to remember, however, that pseu- 
dothermal ghost imaging is typically performed in the 
high-brightness regime. 

It is also instructive to compare our SNR results — 
which apply to reflective ghost imaging through atmo- 
spheric turbulence — to the transmissive ghost imaging, 
no-turbulence results from [5]. The principal differ- 
ences between the SNR behaviors in these two cases is 
as follows. Our SNR expressions have the average in- 
tensity reflection coefficient T(pp) appearing in lieu of 
the magnitude-squared of the field reflection coefficient 
|r(pp)P that is seen in the transmissive case. Moreover, 
the transmissive-case SNR from [5] has neither target- 
induced speckle nor turbulence contributions, so it has 
neither A^Path term nor any turbulence-induced noise- 
magnification factors. The reflective case's A^Source 
term has a factor of (1 -f /3~^) and its A^Mix term 
has a factor of (4/3 + /3^^), but both ji^^ contribu- 
tions are absent from the transmissive case. These ad- 
ditional reflective-case contributions are direct results of 



the rough-surface scattering. In both transmissive and 
reflective operation there is averaging of the source ran- 
domness at the target, but in the reflective case the fields 
scatter incoherently off the target, and their Fourier com- 
ponents are averaged at the bucket detector. This sec- 
ond averaging depends on the ratio of the source size 
to the bucket size, i.e., /3. As /3 increases, the bucket 
detector captures more light, this term becomes less sig- 
nificant, and we approach the transmissive behavior for 
these terms. Finally, the refiective-SNR terms associated 
with the bucket detector measurement in A^Detect and 
A^Mix are scaled by L'^/Ab, relative to the correspond- 
ing terms in the transmissive case. This scale factor is 
the inverse of the solid angle subtended at the target by 
the bucket detector. Its presence in the reflective case is 
a consequence of the quasi-Lambertian scattering from 
the rough surface 

It follows from the SNR differences highlighted above 
that the low-brightness SNR asymptotes for transmis- 
sive and reflective operation only differ by the latter's 
including the solid-angle subtense factor A^/L'^. Their 
high-brightness SNR asymptotes are also similar, with 
the only difference being the reflective case's factor of 
(1+/?^^). The most profound difference, however, is due 
to the A^Path term that is only present for reflective 
operation. This term implies that the SNR of reflective 
ghost imaging has a maximum value, SNRgat , whereas no 
such saturation occurs in transmissive ghost imaging, i.e., 
its SNR can grow without bound as Tj/Tq increases. Fur- 
thermore, in reflective imaging, this limiting effect can be 
quite severe: even with no turbulence, for (3 — 1 we flnd 
SNR,at = 3.26, and for ^ = 2 we have SNR^at = 5.54. 
So, for realistic standoff sensing, the SNR will be limited 
to single-digit values if no further measures are taken 
to suppress the fluctuations arising from target-induced 
speckle and atmospheric turbulence. 



V. COMPUTATIONAL GHOST IMAGER 

We have developed the performance of reflective ghost 
imaging from classical and nonclassical sources in a quan- 
tum Gaussian-state framework. Next we consider a com- 
putational ghost imager, a variant of the pseudothermal 
case in which the "chaotic" signal field is created by de- 
terministically modulating the wavefront of a laser with 
a spatial light modulator (SLM), and no reference-arm 
light beam is required. Instead, knowledge of the SLM 
phase pattern is used to calculate the signal-beam in- 
tensity at the target for use in lieu of a reference-beam 
measurement [7] . If an independent pseudorandom phase 
process is applied to each SLM pixel, then far-field prop- 
agation can be argued to yield Gaussian field statistics 
per the Central Limit Theorem. Thus it is appropriate 
for us to model the computational ghost imager's far- 
fleld statistics with a Gaussian-Schcll model, taking the 
source's coherence length to be approximately the width 
of an SLM pixel, and source's intensity diameter to be 
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the diameter of the SLM. 

The laser that illuminates the SLM will be taken to 
emit coherent-state light, hence the SLM's output beam 
is in classical state and semiclassical theory will give 
quantitatively identical results for its far-field photode- 
tection statistics as the full quantum treatment [21 E [7] . 
Also, because the reference field at the target is com- 
puted, rather than measured, it is most naturally mod- 
eled as a complex-valued Gaussian random process. Con- 
sequently, and without loss of generality, we choose to 
use semiclassical theory for our analysis of computa- 
tional ghost imaging because it lets us deal with only 
classical random processes, rather than a combination of 
quantum field operators and classical random processes. 
We thus use E'g{p,t) to represent the complex envelope 
of the target-plane signal field, and E'j^{p,t) to repre- 
sent the computed complex-envelope of the field that 
would have been present on the CCD camera. These two 
complex- valued, zero-mean, Gaussian random processes 
are completely characterized by their nonzero autocorre- 
lation and cross-correlation functions, which we take to 
be as given in Eqs. ([3|-([5| for their field-operator pseu- 
dothermal counterparts. That these correlation func- 
tions are consistent with classical random-process the- 
ory follows from the pseudothermal source's producing 
classical-state light. 

The principal distinction between computational and 
pseudothermal ghost imaging — which is also its main 
advantage — is that there is neither detection noise nor 
turbulence on its computed reference-arm intensity, 
which is found via 



(56) 



ip{t)^qri / dr Pp{T)h{t - t), 



with Pp{t) — Ap\E'pf{pp,t)\^. The bucket detector pho- 
tocurrent, in semiclassical theory, is given by 



ib{t) 



dr [q7]Pb{T) + Aib{T)]h{t ~ t) 



(57) 



where Pb{t) = J^^dp\E'rp{p,t)\'^ is the impinging photon 
flux, and conditioned on that photon flux the shot noise 
Aib{t) is a zero-mean random process with autocorrela- 
tion function {Aib{ti)Aib{t2)) = q^'nPb{ti)5{ti-t2)- The 



computational ghost image is then obtained from the cor- 
relation function 



/ J-T,/2 



(58) 



Because the shot- noise is zero mean, given Pb{t), it 
does not contribute to {C{pp)), hence the computational 
ghost imager's spatial resolution and image contrast anal- 
ysis is identical to the Sec. |in| treatment of the pseu- 
dothermal ghost imager once we substitute the compu- 
tational system's classical random processes for the pseu- 
dothermal imager's field operators, and then let pB. — > oo, 
and set tr^ = to account for the absence of turbulence 
on the computed reference. Using this approach we im- 
mediately find the spatial resolution to be pL^a. where 



2p| 



2P| 



(59) 



and the image contrast is identical to that given for the 
pseudothermal imager in Eq. ( |32[ ) . In t he l atter regard we 
note that we could replace ip(t) in Eq. (58 1 with the zero- 
mean process, Aip{£) = ip(t) — {ip(t)), thereby achieving 
the same high-contrast operation as is realized with AC- 
coupling in pseudothermal ghost imaging. 

Turning now to the computational ghost imager's 
signal-to-noise ratio, we find more significant differences 
between its analysis and that of the pseudothermal ghost 
imager. In our Sec. |IV| SNR analysis, the second moment 
of the photocurrent cross-correlation function involved 
an eighth-order field moment, two sixth-order field mo- 
ments, a fourth-order field-moment, three fourth-order 
moments of the turbulence, and a fourth-order moment 
of the target's field-reflection coefficient. For the com- 
putational ghost imager's SNR analysis we only have 
to evaluate an eighth-order field moment, a sixth-order 
field moment, two fourth-order moments of the turbu- 
lence, and the fourth-order moment of the target's field- 
reflection coefficient. This simplification is because the 
other terms involved the path and detection noises on the 
reference arm, which are absent from the computational 
configuration. The computational SNR is then given by 



SNR^ 



r'{p,)^e- 



■4(4+0-^) 



A' (1 + , 



^+r^(Pp)^([r- 



r(Pp) 

Ab 



(60) 



It follows that computational ghost imaging shares the imaging, from Eq. (51), with cr|. set to zero because there 
same saturation SNR as pseudothermal and SPDC ghost is no turbulence in the computational imager's reference 
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arm. Likewise, computational ghost imaging has the 
same high-brightness SNR asymptote as pseudothermal 
ghost imaging, from Eq. (52), with cr^ = 0. On the other 
hand, computational ghost imaging has a significantly 
improved low-brightness SNR asymptote, given by 



SNR^, 



(61) 



in comparison with both pseudothermal and SPDC ghost 
imaging. 

For quantitative comparison between the SNRs of the 
pseudothermal, SPDC, and computational systems we 
will first have to put them on equal footing. The dimen- 
sionless brightness quantity I that appears in our SNR 
formulas is the photon flux per spatiotemporal mode of 
the source. However, because of the dramatically differ- 
ent coherence times of the SPDC {^IbTq ^ 1) and the 
pseudothermal and SLM sources {^IbTq ^ 1), compar- 
isons based on equal X values will be for sources with 
dramatically different photon fluxes, whereas we would 
prefer to compare SNRs at equal photon fluxes. To do so 
we introduce the dimensionless quantity Xq = I^bTq ~ 
PpQ/a^fls, which is the photon flux per source spatial 
mode and detector temporal mode. Because all three sys- 
tems are taken to have detectors with identical character- 
istics, identical values for will yield identical photon 
fluxes. Making this substitution we can easily compare 
the low-brightness SNRs from Eqs. ([54|, (55), and (61). 



In the extreme low-brightness regime, the computational 
imager has the best SNR, followed by the SPDC and 
pseudothermal imagers. We can also compare the high- 
brightness limits Eqs. (52) and (53), where the shorter 



coherence time of the SPDC source gives it a clear ad- 
vantage. 

Because the different sources exit their low- and high- 
brightness SNR regimes and enter into SNR saturation 
for different photon fluxes, assessing that behavior re- 
quires us to explicitly choose some operating parameters. 
First, we assume the operating wavelength is A = 1.5 /im, 
the source coherence length is po = 0.15/7rmm, the 
source intensity radius is = 3 cm, and that the bucket 
detector matches the source size, i.e., Ai, = ttaq. These 
give us the on-target parameters of pL — 0.05/7rm, and 
ql = 10 m. The photodetectors are taken to have 77 — 0.9 
quantum efficiency, with the CCD pixel area being Ap — 
0-lpj^- We assume the target is at L = 1km range, that 
it has an effective area of A'j, = 50 m^, and T{pp) = 1 
at the point for which we will perform our SNR evalua- 
tions. All paths are assumed to have moderate-strength 
(C^ = lO^^"' m^/^) turbulence. Finally, we take the pseu- 
dothermal and computational sources to have temporal 
coherence times satisfying Tq = lO^/fis Hz, while for the 
SPDC source Tq = l/lO^l^sHz. 

We consider two scenarios: first, when the source is 
weak enough (Zo = 1) that the SPDC output is still a 
stream of biphotons, but we are not deep into the low- 
brightness regime; and second, when the source is suf- 
ficiently strong {la = 10'') to make the SPDC output 
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FIG. 2: (Color online) SNRs for pseudothermal, SPDC, and 
computational ghost imaging vs. normalized integration time 
Q,bTi. The source brightness is In = 1, so that the SPDC 
output is a stream of biphotons and all three systems are 
well below their respective high-brightness SNR asymptotes. 
Other parameters used in computing these curves are given 
in the text. 
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FIG. 3: (Color online) SNRs for pseudothermal, SPDC, and 
computational ghost imaging vs. normalized integration time 
SlsTj. The source brightness is In = 10'', so that the SPDC 
output appears to be classical and the pseudothermal and 
computational ghost imagers are in their high-brightness lim- 
its wherein they only differ by a factor of e^"^ . Other param- 
eters used in computing these curves are given in the text. 



appear classical . The resulting SNR curves for these 
two cases are shown in Figs. [2]and[3j respectively. From 
Fig. [2] we see that the SPDC imager is the worst per- 
former when it is in its biphoton regime, while the com- 
putational ghost imager is to the top performer there. 
From Fig. [3j however, we see that when the SPDC out- 
put appears to be classical, its much shorter coherence 
time allows it to reach the saturation-SNR limit well be- 
fore the other two systems, whose SNR curves are nearly 
identical. 
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VI. DISCUSSION 

We have derived the key performance characteristics of 
reflective ghost imaging in the context of reahstic standoff 
sensing apphcations. This was done within a Gaussian- 
state framework that aUowed for analysis of classical 
pseudothermal sources and the nonclassical output of an 
SPDC. Our results for the spatial resolution, image con- 
trast, and SNR were analogous to those calculated for 
transmissive ghost imaging in the absence of turbulence 
[21 E] , with a few significant differences that result from 
the roughness of the target surface, turbulence on the op- 
tical paths, and the standoff nature of the measurement. 
First we shall discuss the impact of rough-surface scatter- 
ing in the absence of turbulence, then we shall discuss the 
additional effects incurred when turbulence is present. 

Rough-surfaced targets scatter their impinging illumi- 
nation in a random manner, casting random speckle pat- 
terns in their far fields that, on average, correspond to a 
quasi-Lambertian distribution. Consequently, the ghost- 
imager's bucket detector will, on average, collect only 
a fraction, Ab/L"^, of the target-scattered light. This 
angular subtense factor — which is absent in transmis- 
sive ghost imaging — appears in the average photocur- 
rent cross-correlation functions in Eqs. (22 1 and (23), 



and the noise terms associated with the bucket detec- 
tor, i.e., A^Detect and A^Mix. The surface rough- 
ness has another major effect on the SNR: its time- 
independent speckle pattern implies the existence of a 
finite saturation-SNR value in the limit of long integra- 
tion time. 

The target-speckle contribution that leads to a satu- 
ration SNR appears in A^Path, and is a function of /3, 
the ratio of the source size to the bucket detector size. 
Here we may incur a trade-off between spatial resolution 
and SNR. Specifically, by enlarging the bucket detector 
we can increase the averaging of the target speckle, and 
thus the saturation SNR, without affecting the imager's 
spatial resolution. Ultimately, however, there will be a 
practical limit beyond which the receiving aperture size 
can no longer be increased. Beyond that point /3, and 
thus the saturation SNR, can only be increased by de- 



creasing the source size qq- However, as seen in Eq. (28 1, 



spatial resolution is inversely proportional to oq, so in- 
creasing the saturation SNR in this manner will degrade 
spatial resolution. 

The effects of turbulence on the image contrast and 
spatial resolution of reflective ghost imaging are simi- 
lar those previously found for transmissive ghost imag- 
ing ,llj . In particular: turbulence does not change the 
image contrast; turbulence in the target-to-bucket path 
has no effect on spatial resolution; and turbulence on 
the signal and reference paths degrades spatial resolu- 
tion in the same manner, i.e., degradation occurs when 
they become smaller than the source size. In both SPDC 
and pseudothermal ghost imaging it is therefore advan- 
tageous to propagate the reference field through a con- 
trolled (turbulence- free) environment. Turbulence also 



degrades ghost-image SNR in two ways: first, its logam- 
plitude fluctuations magnify several existing noise terms; 
and second, these fluctuations also contribute to the 
time-independent noise term A^Path, decreasing the sat- 
uration SNR from its speckle-only value. 

In computational ghost imaging, the reference field is 
computed, removing the turbulence on that path, and 
thus improving both the spatial resolution and the SNR. 
The computational case also has significantly reduced 
noise for low-brightness illumination, as there is no longer 
detection noise on the reference arm, making its low- 
brightness SNR behavior superior to that of both pseu- 
dothermal and SPDC ghost imaging. More generally. 
Figs. [2] and [3] show that the computational ghost im- 
ager has a higher SNR than the pseudothermal system 
in both low-brightness and high-brightness operation, 
and than the SPDC system in its low-brightness (bipho- 
ton) regime. However, in high-brightness operation, the 
SPDC ghost imager's significantly shorter coherence time 
enables it to reach its saturation SNR much quicker than 
either the computational or the pseudothermal ghost im- 
agers. Although a cw SPDC will not be able to reach 
that high-brightness regime, pulsed systems are capable 
of doing so [57]. 

Insofar as standoff sensing is concerned, the essential 
conclusion to be drawn from the preceding summary of 
our work is simple: the computational approach is the 
ghost-imaging configuration that shows the most promise 
for this application. It has the best spatial resolution, 
the highest saturation-SNR, and general SNR perfor- 
mance second only to high-brightness SPDC operation. 
In addition, the computational approach obviates the 
need for a physical reference path, which makes it nat- 
urally amenable to obtaining 3D ghost images. In par- 
ticular, pseudothermal and SPDC ghost imaging require 
that reference-arm measurements be made on an inten- 
sity pattern corresponding to the one that is projected 
onto the target. As shown theoretically in [71 and ex- 
perimentally in 8J, computational ghost imaging allows 
reference intensity patterns to be computed at a variety 
of target ranges so that ghost images can be formed for 
these target ranges from the same bucket- detect or data. 
This range sectioning is something that cannot be done, 
without separate measurements for each possible target 
range, in pseudothermal or SPDC ghost imaging. 

We have reserved our final comment for our choice of 
image reconstruction via photocurrent cross-correlation. 
It has already been demonstrated [SI 1^ that ghost imag- 
ing reconstruction can be performed in transmission 
through more advanced methods, such as compressive 
sensing. A preliminary no-turbulence performance com- 
parison between a simple cw-laser radar and ghost imag- 
ing in reflection via cross correlation has shown their sim- 
ilarity in spatial resolution and SNR 12J. If advanced 
image-reconstruction techniques can be successfully ap- 
plied to computational ghost imaging, the resulting sys- 
tem could have notable advantages over current laser 
radars. 
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